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The distribution of the phase angle and the magnitude of the fermion determinant as well as its 
correlations with the baryon number and the chiral condensate are studied for QCD at non zero 
quark chemical potential. Results are derived to one-loop order in chiral perturbation theory. We 
find that the distribution of the phase angle is Gaussian for small chemical potential and a periodic 
Lorentzian when the quark mass is inside the support of the Dirac spectrum. The baryon number 
and chiral condensate are computed as a function of the phase of the fermion determinant and 
we discuss the severe cancellations which occur upon integration over the angle. We compute the 
distribution of the magnitude of the fermion determinant as well as the baryon number and chiral 
condensate at fixed magnitude. 

Finally, we consider QCD in one Euclidean dimension where it is shown analytically, starting 
from the fundamental QCD partition function, that the distribution of the phase of the fermion 
determinant is a periodic Lorentzian when the quark mass is inside the spectral density of the Dirac 
operator. 
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I. INTRODUCTION 



The phase diagram of strongly interacting matter is expected to show several phases as a function of the temperature 
and the baryon chemical potential. Matter in nuclei, in compact stars and in the early universe are in different parts 
of the phase diagram and large experimental and theoretical efforts have been invested to understand their properties. 
Of particular intense interest is the critical end-point. Its existence is expected mainly on the findings of model studies 
that the baryon density is discontinuous as a function of the chemical potential Lattice QCD, which has allowed 
us to determine the nature of the phase transition at zero baryon chemical potential 0], appears to be the natural 
tool to study the non-perturbative phenomena which take place near the endpoint. However, probabilistic lattice 
QCD methods are not directly applicable at nonzero baryon chemical potential: Monte Carlo importance sampling, 
which is at the core of Lattice QCD computations, requires that the Euclidean action is real. At non zero chemical 
potential, though, the quark determinant is complex. This severe obstacle is known as the sign problem. 

Recent numerical progress in understanding the phase diagram of strongly interacting matter at nonzero chemical 
potential has reopened the field. Not only has it been understood that the location of the endpoint in the (/i, T)-plane 
is extremely sensitive to the quark mass 0], it may also be that the dependence of the endpoint on quark mass is 
very different from what was commonly accepted |3|. Because of the sign problem these conclusions where reached 
from analytic continuations of lattice simulations carried out at imaginary values of the chemical potential. Such an 
extrapolation [1,0,H,0| is not without pitfalls. It has recently been demonstrated Q that utmost care should be taken 
when attempting to extract information on the critical endpoint from a Taylor expansion at = 0, [l^ 1^ 3 ■ 



Moreover, it was demonstrated in [T3| that the numerical implementation of the re- weighting approach [isl. IgT 17j is 
extremely delicate even at small values of the chemical potential. 

Lately alternative numerical methods such as the density of states method and the complex Langevin method have 
been explored. Despite early reports of its failure [H, [l^, [2^ , the complex Langevin method has been shown to be 
able to deal with sign problems in simple models and for a gas of relativistic bosons [2l|. On the analytical front, 
the severity of the sign problem was analyzed for QCD at low energy and for models of the QCD partition function 
[13, m, HJ, H^, . The intricate connections between the sign problem, chiral symmetry, and the Dirac spectrum. 



have been understood in the e-regime of QCD [27 1 



In the present work we focus on the density of states method [28|, |29|, |30|, |3l|, |3j, |33[ . In this approach one evaluates 
an observable numerically for a fixed given quantity and thereby obtain the distribution of this observable over the 
fixed quantity. The full expectation value of the observable is then obtained by integration over the fixed quantity. 
This method has had some success when the baryon number, the average plaquette or the phase of the fermion 
determinant is kept fixed. In this paper we are particularly interested in the last approach since it goes back to the 
root of the sign problem. If we would know the exact distribution function of the phase of the fermion determinant as 
well as its correlations with physical observables, the sign problem would have been solved: the delicate cancellations 
due to the fiuctuations of the phase could be realized exactly by an analytical integration over the phase according 
to the distribution function and its correlations. 

We will use chiral perturbation theory to compute the distribution of the phase of the fermion determinant 

/dA|det(D + ^7o +m)|^/e'^/*' e-^™ ^' 

Here 9' refers to the phase of the fermion determinant. It is a function of the gauge field configuration which we 
average over, i.e. exp(2i6'') = dct{D{A) + /X70 + m)/det*(Z)(A) + fijQ + m). Due to the sign problem the distribution 
of the phase is not real and positive. The complex nature, however, is of the simplest possible form: Since 

{S{9-e'))^^=e^'^^f^{6ie-6'))\^^^. (2) 

the ^-distribution factorizes into exp(i6Nf) and a real and positive distribution. Here, Zn^ is the Nf fiavor partition 
function and -^itv/I is the phase quenched Nf flavor partition function. The subscripts Nf and |iV/| refer to averages 
with respect to these two partition functions, in this order. For Nf = 2 this relation reads 

{Si9-9'))^+,=e'">^{S{0-0')h+^., (3) 
^1+1 

where here and below the subscript 1 H- 1 refers to QCD with two ordinary flavors whereas the subscript 1 -|- 1* refers 
to QCD with one ordinary flavor and one conjugate flavor. By definition the fermion determinant of a quark and a 
conjugate quark are each others complex conjugates so that the total measure is real and positive. The ^-distribution 
of the phase quenched theory, {5{9 — 6'')) 1+1*, is necessarily real and positive. Moreover it is normalized to one. Also 
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FIG. 1: A schematic picture of the phase diagram of QCD as a function of the quark chemical potential /x and the temperature 
T. Chiral symmetry is spontaneously broken below the full curve. The dashed curve indicates where the quark mass enters the 
Dirac spectrum. As this happens the nature of the sign problem changes. To the left of the dashed curve the distribution of 
phase of the fermion determinant is a periodic superposition of Gaussians whereas it is a periodic superposition of Lorentzians 
to the right of the dashed curve. We stress that the dashed curve does not indicate a phase transition in QCD. 

the ^-distribution of the full theory {5{0 — is normalized to one. On the r.h.s. of ([3]), however, the ratio 

/2^i+i grows exponentially fast with the volume so that the phase factor, e^*^, must lead to exponentially large 
cancellations. 

As we shall see it is essential to discuss separately the case when is small and the case when the quark 

mass is inside the spectral density of the Dirac operator. We will show below that the real and positive part of the 
^-distribution becomes a periodic superposition of Gaussians when the quark mass m is outside the support of the 
Dirac spectrum. When the quark mass is inside the support of the Dirac spectrum the sign problem becomes much 
more severe [l^, [13 • Figure [T] gives a schematic picture of the phase diagram of QCD as well as the region where 
the quark mass is inside the spectral support of the Dirac operator. As we will show below, the ^-distribution in this 
region is not only very wide/flat, it also changes shape into a periodic superposition of Lorentzians. A hint of this 
dramatic change is already present in ([3]). When the quark mass enters the spectral support of the Dirac operator a 
phase transition occurs in the phase quenched theory while the full theory remains unaltered [13] ■ The exponential 
growth of the ratio Zi^i* / Zi^i with the volume is thus particularly rapid when the quark mass enters the spectral 
support of the Dirac operator. (The link between the Dirac spectrum of the full theory and the phases of the phase 
quenched theory are discussed in detail in [13, [H, [H, [s^] ) . 

The Gaussian shape of the ^-distribution for small ji was first observed numerically by Ejiri in [s^ where it is 
also argued that this form is a natural consequence of the central limit theorem. The change from the Gaussian to 
the Lorentzian form for larger values of /z therefore suggests a breakdown of the conditions for the application of the 
central limit theorem. To cast further light on this we also compute the distribution of the phase, {S{9 — 9')), for lattice 
QCD in one Euclidean dimension. As we will show, it is possible to derive the Lorentzian form of the ^-distribution 
directly from the one dimensional lattice QCD partition function, when the quark mass is inside the support of the 
Dirac spectrum. 

In addition to the distribution of the phase of the fermion determinant we also consider the direct dependence of 
observables O on the phase 6 through the distribution-function 



The integral over 9 obviously gives the full expectation value (O) . The 6'-dependence of the observable shows if severe 
cancellations take place in this integral. Furthermore, the distribution of the observable with the phase allows us to 
address which range of the phase is essential for the full expectation value of O. 

We will compute the distribution of the baryon number operator, its square as well as the distribution of the chiral 
condensate over 9. It is found that the distributions, {O S{9 — 9')), take complex values and that drastic cancellations 



(O 5(9-9')). 



(4) 
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occur when integrating over 6. 

This paper is organized as follows. In section HIl and Hill we briefly recall a few facts about chiral perturbation theory 
which are relevant for the calculation of the average phase factor and the distribution of the phase angle. Then we 
turn to the distributions of the baryon number (section ITV)) . the off-diagonal susceptibility (section |V| and the chiral 
condensate (section IVI[) over the phase angle of the fermion determinant. These one- loop results are all valid for 
/i < m7r/2. Next we discuss the distribution of the phase for an ensemble generated at /i = 0. The difference in the 
phase distribution for e-counting rather than the p-counting in pointed out in section I Villi In section IIXI it is shown 
that the leading order prediction for the ^^-distribution takes a Lorentzian shape for /i > 771^^/2. The Lorentzian form 
is then obtained as an exact result for lattice QCD in one Euclidean dimension in section IXll The remainder of the 
paper discusses the radial distribution of the fermion determinant. 



II. 1-LOOP CHIRAL PERTURBATION THEORY AND THE AVERAGE PHASE FACTOR 

The first step towards obtaining the distribution of the phase is to understand the average of the phase factor. In 
this section we review the calculation of the average phase factor in chiral perturbation theory. 

Chiral perturbation theory [s^ is the low energy effective theory of QCD in the phase where chiral symmetry is 
broken spontaneously. It describes the dynamics of the Goldstone modes, i.e. the pions and the kaons. We shall work 
in the so called p-cxpansion of chiral perturbation theory where the small expansion parameter is 

p^m^^fi^T^^. (5) 

For /i < the chemical potential modifies the pion propagator in the standard way for relativistic bosons. The 

one-loop contribution to the free energy from a pair of charge conjugate pions (the chemical potentials are therefore 
fi and — /i) is thus given by 

Go(m,-m) = -^^logdpL + (p,o - 2*/i)Y), (6) 



where 



Pka = ka integer. (7) 



La 

After a Poisson resummation this can be expressed as [1^ 

Go(m,-m) = -^E/ (0^^^"^""°'° ^og{\f+ml + {po 2zM)'n, (8) 

where the sum is over all integers. The thermodynamic limit is given by the term la = 0. Here, two facts about this 
term^ which we denote by Goly^oo, are essential: i) it is independent of fi ii) it includes the entire 1-loop divergence 
(see [38j for a discussion) . In dimensional regularization it is given by 

Go\v.^ = (-lyni. (9) 

The finite part of the 1-loop free energy, denoted by 5o(m)j contains the sum over the terms with la =/= 0. This results 
in the decomposition 

Go{fi,-fJ.) = Go\v=oo + gn{f^,~fJ')- (10) 

If we wish to keep track of the leading 1/V corrections to the infinite volume result we have to evaluate the sum 
over all four components of the momentum. The finite, fi, L and T dependent part then reads [23| (this expression 
generalizes the result of [i^l for /i = to nonzero chemical potential) 

•^0 a=0 l„ 
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where la runs over all integers and L = (LoLf )^/^. 

When the length of the box is considerably larger than the Compton wavelength of the pion the sum over momenta 
can be replaced by an integral, and the 1-loop contribution to the free energy simplifies to the familiar expression 

.9oM,-M = 2 — V J — cosh — . 12 

n=l 

As the simplest relevant example let us now consider the average phase factor for the phase quenched theory. By 
definition we have that 

(e^^'^'>i+i. = p^. (13) 
^1+1* 

The phase quenched theory in the denominator is identical to QCD at nonzero chemical potential for the third 
component of isospin [4l] |. This has an immediate consequence: since the pions carry isospin charge but no baryon 
charge the free energy of Zi^i* depends on /i while the usual free energy of Zi+i is independent of /i when evaluated 
in chiral perturbation theory. It is this dependence on the chemical potential which makes it possible to compute the 
average phase factor in chiral perturbation theory despite the fact that pions have baryon charge zero. 

For small fi the leading (mean field) term in the chiral Lagrangian, 2m{'ijjip)V , is identical in the two cases and 
hence the phase factor is determined to leading order by the one-loop effect 

(e-')^^^_^g!!:!L=e-(--°)--(-). (14) 

With p-counting (O we have that go{n) — ,9o(a* = 0) ~ Vfi'^T'^ 1 as was discussed in detail in [23 |. 

For fi > a Bose Einstein condensate of pions forms in the phase quenched theory and the mean field terms in 

the chiral Lagrangian contribute to (exp(2i0')). These terms are of order ^^F^V ^ F/L^ ~ L^. Hence, for /i > TO,r/2, 
the strength of the sign problem depends on L even if we scale m,r and /x with L according to p-counting. 

Since the difference of the finite parts of the one-loop free energy appears repeatedly below, it will be convenient 
to introduce the notation 

AGo = AGo(/x, -fJ-, rn, m) = Go(/i, -/i, m, m) ~ Goifi, 'm-i rn) = go{fJ., -fJ., m, m) - go{fJ., /i, m, m). (15) 

Below we will also meet free energies where the chemical potentials are not of opposite sign and where the quark 
masses are different. To be precise we reserve the notation AGq as defined in (fT5)) . and explicitly write the dependence 
on the chemical potentials and quark masses when necessary. 



III. THE DISTRIBUTION OF THE PHASE {fi < 

The distribution of the phase angle can be obtained from the moments of the phase factor [2^ 1 

{S{0-0'))n,^^ e-^P'{e^P'')Mr (16) 



27r 

p—~oo 



The even moments are ratios of a partition function with p additional determinants and inverse conjugate determinants 
and the usual Nf flavor partition function 



(e^-^V. = ^ (J5|§i^^^det-^(Z. + ,,0 + m)\ , 
Zn, XdefiD - fi'jQ + m) / 



(17) 



Since the number of charged Goldstone modes of the partition function in the numerator is p{p + Nf) whereas the 
contributions of the neutral Goldstone bosons from the numerator and the denominator cancel, we obtain 



)n, = e-P^^^+J')'^^^ (18) 



When the quark mass is outside the support of the Dirac spectrum, the contribution to the phase angle of individual 
eigenvalues is in the range [— 7r/2, 7r/2], and we expect half-integer powers of the determinants in ()17p are smoothly 
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connected to results obtained for integer powers. In other words, we expect that the repUca trick [43, l43| can be used 
to analyticaUy continue the moments to half integer values of p. We then find 



P— — 00 



1 o ^^ ?\r2 

2^ 



.iNj0 + \N'lAGa ^ g-m6(-ti2AGo/4 



^g.jv,e+iAr|AGo^3(0/(2^)^ e-'^^o/^)^ (19) 
27r 



After a Poisson resummation this can be rewritten as YM 



{S{0-0'))j, = ^^^e^NfO+lNfAGo y g-(e+2„.)VAGo^ G [-TT, ^] (20) 

VttAGo ^ 



valid for a compact phase angle 9 e [— tt, tt] 
Notice that 



(21) 



so that to be consistent with the general form given in ([3]), the result (|20p shows that the quenched and the phase- 
quenched ^-distributions are identical. Also note that the ^-distribution depends only on AGq. Plots for AGq = 0.2 
and AGq = 10 are shown in figure [21 Notice the different scales in the two plots. For AGq = 10, when the sign 
problem is severe, the normalization to one requires a delicate cancellation. 

As long as the contribution to the phase of the fermion determinant from individual eigenvalue pairs does not 
exceed tt/2 one can unambiguously define the phase of the determinant on [—00,00] as was done by Ejiri (ssj . To 
obtain this distribution simply interpret the angle in (|20p as ranging from —00 to 00. This leads to the Gaussian 
distribution (here for Nj — 2) 

{6(9 0')>i+i = -£^e-o'/^^"+^^", e e hoo, 00]. (22) 

However, when the quark mass is inside the support of the spectrum of the Dirac operator only the phase restricted 
[— TT, tt] can be defined uniquely. We return to this point in section [IXl where we derive the ^-distribution for /.i > m^/2. 

When the angles are noncompact and replica trick can be used, it is useful to represent the ^-function in Eq. (|19p by 
an integral over p instead of a sum over p. Below this will be exploited on several occasions to simplify our expressions. 

IV. THE BARYON NUMBER OPERATOR {/i < m^/2) 

Since the pions have zero baryon charge the baryon number in chiral perturbation theory is automatically zero. We 
will see below that the baryon number at fixed is a total derivative. 

To derive {nBS{0 — 6')) we first compute the correlation between the baryon number and all moments of the phase 
factor 

ins ^ hm ^ C^lfn'' + + ^^^^ ^''^ 

zZi+i d/i det'^(L' — ^70 -|- m) 



To one-loop order in chiral perturbation theory we obtain 

(24) 



\detP{D-f,'yo+m)^^^ \^ + A^TO + HI-)} _ 2p(AGo(-M,A)-AGo(M,A))-p'AGo(-M,/') 



To keep track of the combinatorics it is essential to recall that the one-loop free energy does not depend on the baryon 
chemical potential, that is Go(/^,/i) = Go(/^ = 0). We conclude that 

{UB e2^^^')n-i - - f hm p AGo(-/i,/i)') . (25) 
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FIG. 2: The real part of the distribution of the phase (6(6 - e'))i+i (sohd curve) for AGo = 0.2 left and AGo = 10 right. Also 
shown is the real part of the distribution of the baryon number over 9 (dashed curve). For better comparison the latter has 
been rescaled by (lim/i^^ ^AGo(— /-t, /t)). The fact that the ^-distribution is normalized to unity while the distribution of the 
baryon number over integrates to zero is not easy to see when AGo = 10. This directly illustrates the severity of the sign 
problem. Note that the phase is constrained to ^ £ [— 7r,7r]. 



The delta function 5(6 — 9') is obtained after summing over p. Interpreting the phase angle on (— oo, oo) and proceeding 
in the same way as for the distribution of 9 we obtain 

{UB 6(6-9')),+, = f hm AaGo(-^,^)) (1 + .-^)-^e-^^/^^''+^«°. (26) 

The total baryon number density should vanish because chiral perturbation theory does not include baryonic degrees 
of freedom. This can be seen simply by writing the above expression as a total derivative 

= f hm AaGo(-a., A)) f ^e^^^e-^VAGo+AGo = q. (27) 

The total derivawtive appears because all one-loop contributions to the 2p'th moment of the phase factor are propor- 
tional to p or so that the differentiation to obtain the baryon density leads to an overall factor p. This factor can 
be expressed as a total derivative with respect to 6. Notice that when AGq ^ 1 the extreme tail of the distribution 
over 6 can contribute significantly to the cancellation of the total baryon number. 
If we, as is usually the case, consider the phase on [— tt, tt] we get instead 

{uB 5(6-6')),+, = f lim AaGo(-a.,m)) £ (1 + ^^±^)^^e-(«+-")VAGo+AGo. (33) 

^ ^ n— — 00 ^ ^ 

An illustration of (ns 5(6 — 6')),+, is given in figure [H For small AGq, a small phase angle gives an excess of baryons 
over anti-baryons, which is cancelled by the opposite effect at larger phase angle, resulting in =0. For large 
AGq the plot is quite similar to the ^-distribution which is also shown in this figure. There is however an important 
difference: The integral over 6 of the ^-distribution is unity while the total baryon number is zero. 

The importance of the tail for the cancellation of the total baryon number translates into the importance of the 
terms with large values of \n\. 



V. THE OFF-DIAGONAL SUSCEPTIBILITY {^i < m^/2) 



Even though pions have zero baryon charge chiral perturbation theory gives a nontrivial prediction for the off- 
diagonal quark number susceptibility. To compute this expectation value we start from 

Z,+,(ii,iia) = {dct{D + fi-fa + m)dct(D + fialo + m)). (29) 
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where (. . .) is the quenched average. The average of the ofF-diagonal susceptibihty is then given by 

{x)i+i = i^ 7 r lim (30) 

To one-loop order in chiral perturbation theory we find 

The one loop contribution Go{fia, fJ-b) to the free energy from a charged pion pair made out of quarks with 
chemical potentials fia and Hb only depends on the absolute value of the difference /ia — fib- Moreover, since 
lim^^_,^ d/dfj,Go{fj,, Ha) = we immediately get 

(x)i+i= lim AJ_AGo(^,//,). (32) 

Ma-»M dfl dfXa 



A. The distribution 

To compute the contribution of configurations with a specific phase to the off-diagonal susceptibility we first compute 
the moments (x e^*^* )i+i- We start from 

Zi+i+p|p.(Ata,M&,Ai|Ai) = ( j'^ljpln ^ ™^ det(£' + ^a7o + TO)det(£' + fibjo + m)), (33) 



and evaluate the limit 



(X e^'P" )i+i = lim I .(/Xa,Ai6,Ml/^)- (34) 



For the fermionic Goldstone modes we have an additional minus sign leading to 

Zl + l+p\p-{f^a,f^b,f^\f^) ^ g_pAGo(Ma,-M)-pAGo(Mi.,-M)-p''AGo(M,-M)+pAGo(Ma,M)+pAGo(M6,M)-l-AGo(M„,Mb) 

^1+1 (m,m) 



(35) 



Keeping track of p we find 



(X e2*^'«')i+i = lim 



/(|^AGo(Ma,-M))^ + |^AGo(^.,M) 



g-p(2-l-p)AGo(/J,-/j)^ ^gg-j 



For a non-compact phase angle e oo, oo] we obtain a (5-function in the l.h.s. after integrating over p. Proceeding 
in the same way as for the distribution function of the phase we find 

(X He - ^ ([(1 + i^r + 3^1l|;AG„(.., -rtlL.. + |;5^AG„(..,rt„.,) ^.-.V^-.-.., 

(37) 

The first term between round brackets results from the term ~ p^ in Eq. (|36p which, before summing over p can be 
simply rewritten as second derivative with respect to 0, 

[|;A<^o(... p^^;^e-»-/-..-^ ,38) 

and vanishes upon integration over 9. The 6 dependence of the second term is the same as for the ^-distribution which 
is normalized to 1. Upon integration over the angle 6 we thus recover the expectation value of the susceptibility (|32p . 
Again we emphasize that for AGq 3> 1 contributions from the extreme tail may give important contributions to the 
off-diagonal quark number susceptibility. 
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VI. THE CHIRAL CONDENSATE (^i < m^/2) 



In this section we compute the chiral condensate when the phase angle of the fermion determinant is constrained 
to 6. This quantity is defined as 



(V'V' 5(6 -e')). 



(39) 



Since the chiral condensate is nonzero for T = and /i = this derivation requires also the divergent part of the free 
energy. The required generating functional has different masses 



'^dct''(D-^^■yo+m) 



±gdet(I? + Ai7o + m) 



(det(D + /i7o + ?7i)2) 



(40) 



The desired expectation value (e^'^^ i+i is obtained by taking the derivative with respect to m and subsequently 
the limit m — > m. The combinatorics is much like for the baryon number, but here we have to keep track of both 
mass derivatives and the chemical potentials. This leads to 



1 



1+1 = o ^~ 

I m^ra dm 



(dct(L» + ^7o + m)2) 



(41) 



((V^^Oi+i + ^ - p{Ga{n,~n,rh,m) - Goifi, fj.,m,m)) 
+ (Go(/^, ^J,, rh, rfi) - Go(/x, yU, m, m)) 



-p(2+p)AGo 



Where is the one-loop renormalized chiral condensate at zero temperature and zero chemical potential. For 

p = Q we obtain the one-loop renormalized chiral condensate at nonzero temperature and nonzero chemical potential 



(42) 



which is independent of the chemical potential. 

The distribution of the chiral condensate over the phase 6 is obtained after multiplication by exjp(—ip6) and 
integrating over p 

(-0-0 (5(0 - 6*')) 1+1 = ((V'-0)?+i + hm — {l + ^——){Go{^l,-^l,m,m)-G(i{^,^^,m,rn)) 



m~*m dm L AGo 
-{Goin, /i, rh, m) - Gi^ip., ^, to, to)) 



„2i0 



0FAG^ 



-eVAGo+AGo 



(43) 



The factor 1 + i6/AGo can again be written as a total derivative of the exponential factors. Upon integration the 
contribution from this term vanishes, and we recover the full condensate (1421) 



(V'-0)l + l = 

Again important tail contributions arise for AGq S> 1 



de {iP^ 5{6^e'))i 



(44) 



VII. THE ^-DISTRIBUTION FOR AN ENSEMBLE GENERATED AT /i = 

In the method of Ejiri [l^ one evaluates the ^-distribution as a function of the chemical potential for an ensemble 
generated at zero chemical potential. Here we compute this partially quenched ^-distribution within one-loop chiral 
perturbation theory. 

We start out evaluating the moments of the phase factor for an ensemble generated at zero chemical potential 

{ tTn^'''TW ^^ + - (45) 

^i+i(m = 0) \dct^(L) - /^7o + to) / 
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Wc then obtain the distribution (for /i < 

{5{e~e'))^=o = ^^e'^o. (46) 
VttAG-o 

This one-loop prediction is identical to that for the quenched and phase quenched ensemble: Whether we compute 
the width of the Gaussian for the ^-distribution in the full ensemble generated at ^, or the full ensemble generated at 
/i = 0, or in the quenched ensemble, or the phase quenched ensemble, we find exactly the same result. 

Ejiri also has studied distributions of F = | dct(£' + /i7o + m)|/ det(£'-|-m). His assumption is that the 6 distribution 
remains Gaussian even for a fixed value of F. As we shall see below, this assumption is justified for /i < to 
one-loop order in chiral perturbation theory. 



VIII. THE e-REGIME 

The above analysis suggests that the chemical potential has to be of the order of l/^/V to suppress the correlation 
between the phase and the chiral condensate or baryon density. Such a scaling corresponds to the e-regime [3 
where the dimensionless quantities 

TO EE toEV and ee ^?F'^V, (47) 

are kept fixed for V ^ oo. Here and below S and F are the chiral condensate and the pion decay constant as they 
appear in the chiral Lagrangian. Note that it is possible to go smoothly between the e- and p- regime see [11] . 
In the e-regime the moments of the phase factor remain finite for y ^ oo [l^, US] 

(e^'^^^V, = (1 - 2fL^/rhf(P+^i\ (48) 

where we quote the result valid for m, ft ^ 1 and < to. To obtain the distribution of the phase in the e-regime 
let us rewrite this as 

(e2'P»V/ =e-P^P+^/)'^'^«, (49) 

with 

AGo = - log(l - 2/iVTO). (50) 
It is of exactly the same form as Eq. ((T8)) . So we again find the distribution ((22|) but now with AGo instead of AGq 



{d{0 ~ 0'))l+l = g-9VAGo+AGo_ (5^) 

VttAGo 

The variance of the Gaussian envelope starts out at zero for small fi (i.e. for 2fl^ ^ to) and approaches infinity as 
log(l — x^) for X 1. 
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IX. THE e-DISTRIBUTION FOR > m^/2 

We now turn to the distribution of the phase of the fermion determinant when the quark mass is inside the support 
of the Dirac operator. For low T this means that /i > 771-^/2. 

Because the Dirac operator is only determined up to an operator with determinant equal to unity, 

det[i:» + yU7o + m] = det[A{D + + m)] with det A = 1, (52) 

the sum of the phases of individual eigenvalues of the Dirac operator may differ by multiples or 27r depending on 
the choice of A. When the quark mass is outside the eigenvalue distribution the contribution to the phase of a pair 
X + m, —X + m of eigenvalues is less than tt/2, and in this case the sum of the phases of all the eigenvalues does 
not depend on A. Therefore it makes sense to extend the total phase to (— oo,cx)). For /i > 771tt/2 the phase of the 
determinant, however, differs by multiples of 2tt depending on the choice of A. Therefore, when /i > m-^12^ it only 
makes sense to define the phase modulo 27r. 

As before, the (5-function, 5{6 — 9'), will be obtained from the moments of the phase factor which now are dominated 
by the leading order term in the chiral expansion. Not surprisingly, this leads to a much wider ^-distribution. What is 
perhaps somewhat surprising is that, as will be shown below, the distribution now takes a Lorentzian shape. Because 
of ambiguities in the phase angle we do not expect that we can use the replica trick to calculate half-integer moments 
of the phase factor. Therefore we will only evaluate the even moments, {exp(2ip9')) , with integer values of p. This is 
sufficient to obtain the full distribution of the total phase angle, 29 € [— tt, tt], of det(L' -I- fijo + m)^ relevant for the 
two flavor theory 

(,5(26' - 26*')) = - y e-2'P«(e2*P«'). (53) 

P— — 00 

Alternatively this can be seen as the combination ^[{S{9 — 9')) + {6{9 — 9' + 7r))] of the distribution of the phase angle, 
9, of det(D + /i7o + m). 



A. Bosonic partition function 

The moments of the phase factor involve inverse powers of determinants, c.f. Eq. p7p . As was realized when 
investigating the partition function with one bosonic flavor such inverse determinants lead to a phase transition at 
fi = m7r/2. In order to compute the moments of the phase factor for fi > m.,^12 to leading order in chiral perturbation 
theory we therefore first recall the explanation of the exact results for the bosonic partition function (obtained 
by inte grat ion over the Goldstone manifold ^\ or from the Cauchy transform of the fermionic partition function 
[isl-liollsot) in terms of a mean field argument. 

The observation of [i^ is that the bosonic partition function 



det(Z? — //70 -I- to) 



^ det(£' -f /x7o to) / \ det(L' + ^70 -|- to) (D - /Z70 -f to) 
at a mean field level behaves like 

(dct(£' - ^70 -f to)) 
(det(L' + ^70 +rn){p - /X70 + to)) ' 



(54) 



(55) 



The reason is loosely speaking that the inverse determinant must be regularized in order to be convergent and that 
Grassmannian mean field terms are absent. 

The denominator of Eq. (|55[) is the phase quenched theory which has a phase transition at /i = to^/2. The mean 
field result for the phase quenched theory is given by 

(det(L> /^7o + m){D - /Z70 + to)) = e"^^^ (56) 

where [H 
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is the static Lagrangian for /i > 771^^/2. The average of the determinant in the numerator of Eq. (j55p is the famihar 
one flavor partition function (which is independent of /i in chiral perturbation theory) 

(det(i:' - /i7o + m)) = e-^^°/2, (58) 

where 

Lo = -2mS (59) 

is the mean field Lagrangian at /.( = 0. In conclusion, the mean field result for the bosonic partition function is given 
by 

{—-—A -) = ^-VL,/2+VL,_ (go) 

Met(D + ^7o + m)' ^ ' 

As shown in detail in [13, IS] this gives the correct mean field physics. Note the striking difference with the fermionic 
partition function (j58p which is independent of the chemical potential. 

B. The quenched ^-distribution 

Let us now use what we learned from the bosonic case to compute the quenched distribution of the phase of the 
fermion determinant for > We will first show that 

(e2»ps') = e-^^«lPl, (61) 

where Lb = Lq — Li with Lq and Lj given above (note that Lb > 0). 

Since by charge conjugation symmetry (exp(2ip0')) = (exp(— 2zp0')) this expectation value only depends on the 
absolute value of p, and we only need to consider p > 0. First, we rewrite the moments as 



,2ipe'\ _ I (det(£' + [i-jQ + m) det(j:' + /i7o + to))p 
' \ {Ac\.(p + + to) det(£) - ^70 + m))P 



(62) 



Now the contribution from the denominator is the inverse of the replicated phase quenched theory. This was worked 
out in 

\ = eF^^^ (63) 

((det(i:> + /i7o + m)det(i:i - ^70 + to))p) 

The contribution from the numerator is just 

((det(i:' + ^70 + m) d,cl{D + m7o + to))^) = 6"^^^°, (64) 

which together with the previous result reproduces Eq. (j6ip . The sum over p results in 

(^(26* - 261')) = - V e-2''''Pe-^^«IPl 

P— — OQ 

= i V 2^^^ (65) 

Tx {VLbY + {29 + 27rn)2 ' ^ ^ 



The sum over n can be evaluated as 



This is a compactified Lorentzian, centered at zero. We recall that 20 E [— tt, tt] is the phase of det(I? + /i7o + to) 



2 
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C. The unquenched ^-distribution 

To calculate the unquenched ^-distribution function we again consider the moments {exp{2ip9')) Nf ■ They can be 
rewritten as 

{e^'^'^N, = ^( /'nfnY^° \"'\ + + + + (67) 

Znj- det (D + ^jo + fn) 

where we have introduced u ~ p + Nf /2. By charge conjugation, this expectation values does not depend on the sign 
of u = p + Nf/2, i.e. it only depends on \u\, and it only has to be calculated for p > —Nf/2. We separately consider 
the cases p> and 0> p> —Nf/2. 

For p > there are inverse powers of dot*, and we apply the rules of section flX Al 



^(e2'f^'det^^p + M7o + m)) = ^( 



+ /.t7o + m) \ 



{dci{D + /i7o + m) dct(D - ^70 + m))P 

1 (det'P+^^p + /^7o + »^)) 

~ {{dci{D + /i7o + m) Aet{D - ^70 + m))p) ' 

where the final equality holds at the mean field level. The contribution from the denominator follows again from the 
result of the replicated fcrmionic theory ((dot dct*)P). see Eq. (|64l) . The numerator is equal to cxp(— (p + Nf/2) Lq) 
and the normalization, 1/Zmj, gives cxp{Nf/2 Lq). Therefore the Nf dependence cancels, and wc find the quenched 
result for p > 

{e^ve')^^ ^^-pVLs^ p>0. (69) 

Here, we extended the equality to p = which is satisfied trivially. 

Now, let us look at negative values of p. This means that det* is in the numerator and det is in the denominator. 
For —Nf/2 < p <Q the moments can be rewritten as 

1 dct^ {D + fi-fo + m) ^^^N 



Zn; det* ^(D + ^70 + m) 

'-dei"'[D + ii-io + m)), 



1 det* IPl(L> + ^7o + m) ,_^jv, 



Zn/ det'*'' {D + //70 + m) 
= -^{{dct{D + ^i^o + rn)dct*\D + ^l^^i+m))\P\dci^^-^\^\{D + ^l-|o+m)). (70) 

Note that both exponents are positive. The \p\ pairs of conjugate quarks form a pion condensate while, at the mean 
field level, the Nf — \p\ quarks are passive spectators resulting in the average phase factor 

(e^'P'^Nj = e^T^io-vHL,-y(7V,/2-H)Lo^ -Nf/2<p<0 

— gV\p\{Lo-Li) ^V\p\Lb ^ j-j-j^-j 



Note that it smoothly connects to the p > result ([5^ . 
Combining the above results wc find 



(72) 



for any integer values of p. 

The general result (jTS]) implies that the ^-distribution is a Lorentzian times the phase factor: 



{S{2e - 20')>i+i = - y ^~2^ep^-VLsi\p+l\-l) ^ e^'^ ^"""^" wf."^^^^^''^o.^ ■ (73) 
^ ^ -n ^ TT cosh(V^LB) - cos(26') ^ ' 



1 ^ „~2iep„-VLB(\v+\\-l^^„1^e^ sinh(yiB) 

p— — 00 

As wc have seen previously, the unquenched 1+1 distribution is related to the phase quenched 1+1* distribution 



(<5(20 - 20'))i+i = e^^'^^me - 20'))i+i*- (74) 
^1+1 
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Comparing this with Eq. ([55)) and Eq. (|73p we sec that the quenched and phase quenched 6'-distributions are identical 
also for fi > 

In conclusion, we have shown that the ^-distribution is nonanalytic at the point where the quark mass enters the 
support of the Dirac spectrum. This implies, for example, that the distribution of the phase in this regime cannot be 
obtained by analytic continuation from imaginary values of fj. (see j53l [53 | for a discussion of the analytic continuation 
of the phase factor to imaginary values of the chemical potential). 



X. THE DISTRIBUTION OF THE BARYON NUMBER AND THE CHIRAL CONDENSATE {n > m^/2) 

In this section we compute the distribution of the baryon number and the chiral condensate over the phase angle. 
As for the distribution of the angle itself we will work to leading order in chiral perturbation theory which is the mean 
field result for /.t > 



A. The baryon number 

In order to work out {nBS{2d — 26')) we need the moments 

■dctiD + jl-fo + myf), (75) 



1 ,detP(i:' + ^t7o + m) 



Z^f det^ (Z? - ^7o + m) 

where the chemical potential for the Nf quarks is denoted by //. The distribution is then obtained after differentiation 
w.r.t. p. at fl = fi, multiplication by exp(— 2ip0) and summation over p. 

For p > the bosonic mean field rules discussed in previous section lead to a factorization of the moments as follows 

1 (det'P(i? + /^7o + m)det(D + /i7o + m)^/) 

(76) 



— ^7o + m) det^{D + ^70 + to)) 

Since |/i — /i| < to^ there is no condensation of pions for the partition function in the numerator. It follows that there 
is no dependence on fl and hence all terms with p > vanish after differentiation w.r.t. p.. 

When p is negative the det'''' {D — /i7o + m) is in the numerator and condensation of pions occurs. This leads to a 
dependence on [1 through the mean field Lagrangian 

oy2 2 

Lji-t,, m) = -2F2(^ + ;i)V4 - ^^^^^2p2 - (77) 

Note that this reduces to Lj given in ([57]) for p. = fJ.. 

As in the previous section we must consider separately the cases —Nf/2 < p < and p < —Nf/2. For —Nf/2 < 
p < the moments are given by 

— {e^P^'>'(^^)det''f{D + fljo + m)) = e~^\p\VL,{-^...^.)+[p\VL,i-^.,^.)+\p\VLo (^g) 

at mean field level. While for p < —Nf/2 we find 

^(e2p.e'(M)^^^iv,(^^ ^g-jv,yL,(-^,A)+b|VL.(-M,M)-(bl-JV,)VLo. (79) 

In both cases the derivative w.r.t. p pulls down the prefactor V[d/dil]Li{—fi, jl) but multiplied with a different 
numerical factor. This leads to 



1 



TT 



d 



nB5{2e - 29'))Nf = - -p:Li{-ii, fl) (80) 
dp 



^ (_2)|p|e"^'P^e""^^^('P+^^/^'"^^/^) + ^ (^_]^^-^^-2ip9^-VLB{\p+Nf/2\-Nf/2) 
-Nf/2<p<0 p<-Nf/2 
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For Nf = 2 there is only one term in the first sum, namely p = — 1, and it can be included in the second sum 



(715,5(20 -20')) 1+1 = - 



VL, 



(-2) E ^ 

p<-i 



-2ipe ^~vlb{\p+i\-i) 



The sum can be performed analytically, 



(nB<5(20-20'))i+i 



VLi 



,VLe 



2 rVLi 



32^^- — — log(e^^^ 
2i dd 



(81) 



0- (82) 



The total baryon number density is given by the integral over the distribution (recall that 29 S [— tt, tt]) and vanishes. 
The distribution of the baryon number over the phase angle is proportional to a total derivative but not of the 
distribution of the phase as was the case for fj, < rnT,/2. 



B. The chiral condensate 

As in section IVTl we now dentoe the mass of the Nf quarks by m and differentiate with respect to this mass. The 
computation is somewhat analogous to the one given in the previous section except that the terms with positive p 
also contribute. For Nf = 2 we find 



(V^V 5{29 - 26')),+, = t±^y ^-2:pe^-VLs{\p+i\-i) _ ± _^l,(^, ^) ^ g-2.p9g-yL,(|p+i|-i) (33) 



d 



din 



-1 



p—O P— — 00 

where 



8Ai2F2 

The sums can be rewritten as 



LJ{m,^)^^2^.'F'^^Li^^±fl. (84) 



2SF e-^^^ 2Y.V 2 Y?m 



e 



The first and the last term both vanish upon integration over 9 e [— 7r/2,7r/2] and leaves, 2FE, which, after dividing 
by the volume, is the expected mean field value of the chiral condensate for Nf ~ 2. Note that the amplitude of the 
first term is exponentially small while that of the last term is exponentially big. The severe cancellations which take 
place upon integration of the last term over 9 are just like those for the baryon number. 



XI. QCD IN ONE EUCLIDEAN DIMENSION 



In this section we will show that for one dimensional QCD the distribution of the phase of the fermion determinant 
changes from Gaussian to Lorentzian shape when the quark mass enters the Dirac spectrum. 

Lattice QCD in one Euclidean dimension (time only) with gauge group U{Nc) is sufficiently simple that we can 
solve the partition function and moments of the phase factor analytically starting from the fundamental partition 
function. The reason is twofold. First, there is no Yang-Mills action, and second, the staggered Dirac operator, M, 
can be reduced to the determinant a N^ x N^ matrix [55| 

detAf = 2-"^'^ det[e"^^ + e""^"^ + e'^fU + e-"^C/'f], (86) 

where U G U{Nc). The analogue of m^/2 or rriN/Nc is fic ~ sinh^^ m and n is the number of lattice points. The 
eigenvalues of M are located on an ellipse of width sinh/^ along the real axis. This means that the quark mass is inside 
the eigenvalue domain when fi > ^c- In the limit nNc 00 the ratio of the full partition function and the phase 
quenched partition function approaches one for the SU{Nc) theory whereas this ratio the U{NA partition functions 
shows a phase transition when the quarks mass enters the Dirac spectrum exactly as in QCD j56| . For this reason we 
study the U (Nc) lattice model rather than the SU{Nc) lattice model. In adddtion, the U (Nc) model is mathematically 
simpler than the SU{Nc) lattice model. The partition function is defined by 



ZNfificfJ.) ^ / d?7detM. 



(87) 
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Despite its simplicity many interesting things can be learned from QCD in one dimension. For example, in [56| it was 
found that spectral density of the Dirac operator is a highly oscillatory function when the quark mass is inside the 
ellipse of eigenvalues, and that the link between these oscillations and the chiral condensate is exactly the same as 
what was found for 4d QCD with dynamical quarks [27{ . 

Below we will show that the distribution of the phase of the fermion determinant in one dimensional QCD also 
undergoes a transition from a Gaussian to a Lorentzian shape when the quark mass enters the eigenvalue spectrum. 
For simplicity we only work out the quenched distribution. As above we start from the moments of the phase factor 

Notice that the expectation value only depends on \p\. Following [5^, where the first moment {p = 1) was worked 
out, we rewrite the ?7-integral as 

2.pe'. ^ f detP(l - £/e»A^-"^OdetP(l - [/te-»M-nMc) 

^ ^ 7y(jv<,) detP(l-i7e-"A'-"fc)detP(l-t/te"A'-"^''^)' ^ ' 

In this form the Conrcy- Farmer- Zirnbauer formula [57| can be applied directly for /i < /ic- In the large limit the 
result simplifies to 

2 



(^^^^p, ) = {e''<^ )P = (^1 _ £_ j . (90) 

If /ic is interpreted as the chemical potential for which the quark mass enters the eigenvalue domain, this is exactly 
the same form as we obtained for the e-regime of QCD in section I Villi Hence we find the expected Gaussian form 
for the distribution of the phase of the fermion determinant, 

{d{9 - 9')) = -^e"^'/" for /x < /ic, Nc ^ oo, (91) 

where n = ~ log(l — /i^//i^). 

For fi > fj,c the conditions for applying the Conrey- Farmer- Zirnbauer formula directly are violated. Now, however, 
we instead can rewrite the determinants containing as 



det^'(l - [/te-"A'-"M<=) ^ ^_2pnN,f, detP(l - [/e"A'+"Mc) 



det^'(l - [/te"^'-"^0 det^'(l - [/e-"M+nMc) ' 



(92) 



so that the entire integrand in Eq. (|89p only depends on U . This implies that when we expand the denominator in U 
(which is allowed for fj. > fie) only the constant term is nonzero upon integration over U . Using that the moments of 
the phase factor only depend on \p\ we obtain the exact result 

(e2»P»') = e-2"IPl^=^, for fi > fi,. (93) 
Summing over p in Eq. (|53p results in the compact Lorentzian c.f. Eq. (|65p 

{6(26 - 29')) = - , ^^M'^^N.fi) fi> fx, , 29eh7:,n]. (94) 

^ ^ ^' TT cosh{2nNcn) - cos(2e') t-c , l , j \ j 

We stress that this exact result is valid for any value of Nc- 

Note that we have computed the distribution of the phase angle of the square of the fermion determinant, i.e. of 29. 
The reason is that this does not require the use of the replica trick. By comparing the numerical result for half-integer 
moments with the analytical result ((M]) obtained for integer moments one finds that the replica trick does not work 
when quark mass is inside the eigenvalues. See figure [3] The only exception is the case /ic = 0: Then the rewriting in 
Eq. ((92|) results in 2p powers of the determinants |58| which are then well-defined for half-integer p. The expression 
for the odd moment when /ic = is therefore also given by ()93|1 . 
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FIG. 3: Numerical evaluation of the quenched moments of the average phase factor in one dimensional QCD versus p for 
/ic = 0.1, n — 4 and Na = 3. As indicated by the lines the even moments join smoothly in accordance with (|93p . However, the 
even and odd moments are not smoothly connected for < /Xc < M- 



XII. DISTRIBUTION OF / = log | det(I> + ^70 + m)|/det(I> + m) FOR fi < m^/2 



So far we have considered distributions of the phase of the fermion determinant. As we now show it is also possible 
to compute the distributions as a function of the absolute value of the fermion determinant. We will do this to one 
loop order in chiral perturbation theory using the replica trick. Since only the case fi < will be considered there 

are no issues with the use of the replica trick. 

Since the absolute value of the fermion determinant depends on the large eigenvalues of the Dirac operator we 
analyze the distribution of / = log[| det{D + fi^Q + m)|/ det(_D + m)] which, as wc shall see below, depends only on 
the finite difference of the one-loop free energy at fi and at = 0. In [s^] the distribution of F = exp(/) was studied 
in lattice QCD using the Taylor expansion method. Since determinants fluctuate by many orders of magnitude we 
feel that it is more appropriate to analyze the distribution of the logarithm of their magnitude instead. The two 
distributions are related by a simple transformation 



{S{f-n)=F{S{F-F')), 



(95) 



where /' is the magnitude of the logarithm of the ratio of the determinants - its fluctuations are induced by the gauge 
field fluctuations. 

To compute the distribution of the magnitude of the logarithm of the determinants we rewrite the ^-function as 



m-n 



1+1 — 



^/g-'p(/-/')\ 



dp _i f ( det(£) + /i7o + m) dei{D ~ /i7o + m) \ ^ ..p, , , A .or\ 
\\ det\D^m) ) det(I? + ^7o+r«)^. (96) 

For even ip the average can be interpreted as a partition function with bosonic and fcrmionic flavors. We will 
calculate this partition function to one loop order in chiral perturbation theory. Since we consider the magnitude of 
the determinant we expect that the moments will be analytic in p and can be analytically continued to imaginary ip 
As far as we know this is the first case where the replica trick is used this way. 

Using the same one-loop combinatorics as before, we find after analytical continuation to imaginary ip, 



m /')) 



1-1-1 



2TT 



(97) 



where 



£;/ =2AGo(m)-4AGo(m/2), 



(98) 
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FIG. 4: The /-distribution in 1-Ioop chiral perturbation theory. The distribution of the partition function with 
/ = log(| det{D + fi'yo + m)]/ det(_D + m)) in a box with Vm'^ = 10. The temperature is fixed and as /i increases the distribution 
becomes broader and moves away from zero. 



= iAGo(M)-2AGoW2). 
The integral over p is Gaussian and can be evaluated by completing squares. This results in 

1 ^ 



cr/ V27r 



(99) 



(100) 



Both Ef and are positive. In the thermodynamic limit at nonzero T and fi we can see this using Eq. (|12p 



4 



7^2 „2 



71=1 



,2 fin , /M". 

cosh( — )-4cosh( — , 



Tr2 Z-^ r)2 ^ OT^ / 



27r2 



'2T' 



(101) 



Similarly we can write Ef as 



Ef = 



7r2 2^ 



cosh(f^) 



2cosh(^) 



2 — ^5 — TT — cosh( — )(cosh( — ) - 1) 

Tl=l 



T 



T 



(102) 



Exactly the same combinatorics can be applied to the finite L expressions for and Ef (sec (jlip ) resulting in the 
positivity of these quantities at finite L and Lq. 

Let us make a simple cross check of the formula for Ef and dy. Since / = for = the expectation value and 
the variance of the /-distribution must vanish in the limit fi which is indeed the case (see figure |4]). 

In oder to better understand the structure of the result it is useful to work out the combinatorics for an arbitrary 
number of flavors Nf 



1 



U-NjEj/2) 



(T/V StT 

We note that, an increasing number of flavors simply shifts the average value of /. 



(103) 
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A. The distribution of the baryon number over / 

Even though the baryon number is zero when evahiated in chiral perturbation theory it does not necessarily vanish 
when evaluated for a constrained fermion determinant. In scction lTVl wc derived the distribution of the baryon number 
for fixed phase. Here we compute the distribution of the baryon number as a function of /. 

In order to compute {nBS{f — we denote the chemical potential in the usual two flavor determinant by fl 

instead of /x, then differentiate with respect to fl and finally take the limit jl ^ fi. The 5 function is represented as in 
the previous sections 

(ns <5(/-/')>i+i (104) 

= ^ hm ± r ^e--/ / fdet(D + M7o+Hdet(D-,,o + ^)\ - ^ \ 

2Zi+i A-P d/i 27r \V <let\D + m) J V ^ /u y 

To one loop order in chiral perturbation theory this becomes 

{nB 5{f-n)i+i = [;^(Go(-^,m)-2Go(0,A))]. , T ^ e-^'(/-^/)-^'^>^ (105) 



The Gaussian integral over p results in 



{ns Sif /')>i+i = -[^ (Go(-/i,/i) - 2Go(0, A))]^^^^^^(,5(/ - /'))i+i. (106) 

The baryon number operator is not positive definite and neither is its distribution over /. It changes sign at the 
expectation value of the Gaussian distribution so that the total baryon density vanishes 

POO 

(715)1+1 = / df {uB S{f - /'))i+i = 0. (107) 



As is the case for the distribution of the baryon number over 9 the zero value can also be obtained by noting that the 
integrand is a total derivative. 



B. The distribution of the chiral condensate over / 

In this section we derive the distribution of the chiral condensate over /. As above we represent 5{f — /') by an 
integral over the moments so that 

{^P^P ^(/-/')>i+i 

= ^ lim A r d_P,r^rf I /det(i^ + ,7o + „.)det(7^-M7o+H\ ^ .^^.^^ ^ ^ ^-„A (^^3) 
The combinatorics of possible one-loop contributions of Goldstone bosons leads to 



(ViV"5(/-/'))i+i (109) 

dp r d 
— \ ip-r-:r 
27r dm 



[*PT^ iGo{^J■,^J','m',7h) + Go(-Ai,^,TO,m) - 2Go(0, ^i, m, m)) 



+4^ {Go{0,m) - Go{0,m))]^^^^^^e'-^-y+^P(^f-f\ 
= ('0V')i+i('5(/ - /')>i+i - ^[Go(Ai,/x,m,m) + Ga{-n, fi,m,m) - 2Go{0, fj.,m,m)]^^_^^-^ — {S{f - ./')>i+i- 

The first term in the last line gives the chiral condensate upon integration over / while the second term in the final 
line integrates to zero in precisely the same way as in the case of the baryon density. 
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C. The /-distribution evaluated in an ensemble generated at ^ = 

In [3^ the distribution of = | det(D + /.J7o + m) |/ det(Z3 + m) is studied in lattice QCD for an ensemble generated 
at ^ = 0. We will again study the distribution of / = logi^ for this case. It is given by 

^ {Sif-r)dct\D + m)) (110) 



Z,+i{fi 0) 

= / , f ^e-'^/ / f ^^^(^ + + ^^^^^ -,io+m)Y ,,,.(^ ^ 

^i+i(M = 0)i_oo2^ \V det'(i? + m) J ^ 7 

When evaluated to one loop order in chiral perturbation theory we find 

-^{S{f - /')dct^(i? + m)) ^ r ^e-'P^f-'^f^-^^'fP' = -J—e-^f-'^fy'/i^^^l (111) 
Zi+i J-00 27r afV^n 

where 

Ef = 2AGo{pL/2). (112) 

In comparison to (|100p we see that only the expectation value of / has changed whereas the variance takes the same 
value as in previous sections. 



XIII. CONSTRAINING BOTH 6 AND F FOR (fi < m^/2) 

In oder to understand what happens if both the phase and the magnitude of the fcrmion determinant arc fixed we 
need to compute the correlation between the moments of the phase factor and F. 
Let us consider the correlation of any moment of the phase factor and F 

I A&i^{D + /x7o + m) det''(D + /i7o + to) det''(£) - /^7o + to) 
\det^(D-^7o + m) A<A'^{D + m) det«(D + to) 

/ det*'(Z? + /i7o + to) \ / dct'^(Z? + /170 + to) dct'(D - /i7o + to) 



\ defP (D - //70 + to) / \ det'' (L» + to) det'' (£> + to) 

^ g-p=AGo+9^AGo-VAGo(A'/2) _ g-p2AGog9^AGo-4g^AGo(A'/2) ^ q ^^^2.3) 

The reason is that terms linear in p in the first exponent cancel completely. In other words, even though there are 
bound states (Goldstone bosons) with non zero charge which potentially can couple the phase factor to the absolute 
value of the determinant, their contributions exactly cancel each other. This is also the case if the average is calculated 
for N j dynamical flavors. 

We have thus shown there are no correlations between the absolute value of the fermion determinant and the phase 
to one- loop order in chiral perturbation theory (for < to^/2). Hence we automatically find 

m - mo - o'))i+i - {s{f nh+i{s{o ~ e')),+, (im) 

and 

{nB5{f - f')5{6 - 0'))i+i = {nB5{f - f))i+i{5{6 - e'))i+i + {5{f - f )),+i{nBS{9 - 0'))i+i. (115) 
One can convince oneself that this factorization does not hold for /.j > to^/2. 



XIV. CONCLUSIONS 



The distribution of the phase of the fcrmion determinant for QCD with nonzero quark chemical potential has been 
computed to leading order in chiral perturbation theory. When the quark mass is outside the support of the Dirac 
spectrum (small ^) the distribution becomes Gaussian whereas the distribution is Lorentzian (modulo 2tt) when the 
quark mass is inside the support. This non- analytic behavior is also found for QCD in one Euclidean dimension by a 
direct evaluation of the involved partition functions. 
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The distribution of the baryon number and the chiral condensate as a function over the phase angle has also been 
computed in chiral perturbation theory. The results show analytically that extreme cancellations are essential for the 
vacuum expectation values of these fundamental quantities. 

The ratio of the magnitude of the fermion determinant to its value at = is ultraviolet finite and can be studied 
within chiral perturbation theory. We have computed the distribution of the logarithm of this ratio, /, as well as the 
distribution of the baryon number and the chiral condensate over /. Contrary to the ^-distribution the distribution 
of / is real and positive. In fact, within one-loop chiral perturbation theory for /i < there are no correlations 

between the phase and the absolute value of the fermion determinant. 

The results obtained here are complementary to lattice results obtained by Ejiri [3^. The results for one-loop chiral 
perturbation theory when the quark mass is outside the eigenvalue distribution of the Dirac operator, confirms the 
Gaussian shape of the ^-distribution first found in lattice simulations (ssj . The analytical results, however, also show 
that exponentially large cancellations may take place when integrating over d. Not only are these cancellations essential 
in order to measure the baryon number and the chiral condensate correctly, the extreme tail of the distribution may 
contribute significantly to the final result. A small non Gaussian term in the tail of the ^-distribution therefore could 
be the dominant term after integration over d. The precise form of this tail is of course difficult to access numerically. 

The Lorentzian shape of the distribution of the phase valid for larger values of the chemical potentially shows that 
one should not take for granted that the conditions for the central limit theorem are satisfied. The nonanalyticity 
means that the Lorentzian shape cannot be obtained by analytic continuation from imaginary values of the chemical 
potential. Since the Lorentzian form is present also for quenched QCD this prediction can be tested in lattice QCD 
without worrying about the sign problem. Numerical convergence is however expected to slow because of the large 
ffuctuations of the phase. If staggered fermions are used one also has to address the issues raised in (soj . 

Finally let us stress that both the Gaussian and the Lorentzian forms for the ^-distributions found here are leading 
order predictions of chiral perturbation theory. It would be of considerable interest to work out the next to leading 
order corrections. It seems natural that these terms will give corrections to both the shape of the ^-distribution and 
to its width. 

The analytical work of this paper was inspired by new developments in the numerical density of states method. Such 
interplay between numerical lattice QCD and analytical methods, is essential for progress towards our understanding 
of strongly interacting matter. In this paper this was illustrated by simulations of one-dimensional lattice QCD. Even 
if the analytical results do not yet offer a direct solution of the sign problem, they allow us to better understand the 
regions where current numerical methods can be applied [b^]- 
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